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Abstract 

In the MHD description of plasma phenomena the concept of magnetic hehcity turns out to 
be very useful. We present here an example of introducing Euler potentials into a topological 
MHD soliton which has non-trivial helicity. The MHD soliton solution (Kamchatnov, 1982) 
is based on the Hopf invariant of the mapping of a 3D sphere into a 2D sphere; it can have 
arbitrary helicity depending on control parameters. It is shown how to define Euler potentials 
globally. The singular curve of the Euler potential plays the key role in computing helicity. 
With the introduction of Euler potentials, the helicity can be calculated as an integral over the 
surface bounded by this singular curve. A special programme for visualization is worked out. 
Helicity coordinates are introduced which can be useful for numerical simulations where helicity 
control is needed. 



1 Introduction 

Magnetic helicity is a topological characteristic of magnetic field structures which includes the twist- 

■ 

'.^ ing and the kinking of a flux tube as well as the linkage between different flux tubes (Moffatt, 1978, 
^ ; Biskamp, 1993). Among its numerous applications are dynamo theory (Moffatt, 1978), investigation 
of magnetic reconnection (Wiegelmann and Biichner, 2001), theory of relaxation (Taylor, 2000), and 
even the collimation mechanism of astronomical jets (Yoshizawa et al., 2000). Magnetic helicity is 
defined volume integral 

K = [ (A-H)rfy (1.1) 
Jn 

where B is the magnetic field and A is the vector potential 
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B = VxA. (1.2) 

Helicity (|1 . 1| ) is gauge invariant, because under the transformation A' A + V0, it then is 
changed by 



5K = [ i\/(p-B)d^x = (f (p{B- dS) =0, (1.3) 
Jq Jdn 

if Bn\dn = 0, where n is the vector normal to the boundary dQ. For Bn\dQ 7^ the surface 
integral does not vanish and the helicity becomes gauge dependent. Generally speaking, there is 
the possibility to define the helicity for difference between the original field and the vacuum field 
(Schindler et al., 1988; Biskamp, 1993; Priest and Forbes, 2000) which helps to give the helicity a 
physical meaning for more realistic conditions. Nevertheless, we will restrict our consideration to the 
classical case -Bnlen = 0, leaving a more general definition of the magnetic helicity for future studies. 
When the Euler potentials a, P are used. 



B = Va X V/3, (1.4) 
there is the following problem related to helicity. It can be easily verified that 



A = -pVa (1.5) 

(or A = aVP) is the vector potential ( |1.2D for the magnetic field ( [1.41 ). Then helicity vanishes at 
the level of the scalar product (A • B) = 0. It is known (see, for example, Biskamp, 1993) that the 
vector potential can be presented in the following form (Clebsch representation) 



A = -pVa + Vij, (1.6) 

where the function ip (contrary to in (|1.3|) ) must be multi- valued. This implies that the function 
ip has a surface Sj inside the volume Q where it has a jump, then the contribution from the jump 
surface Sj is added to the integral over dQ in equation (|1.3| ) which results in the nonzero helicity. 
The solution to the questions how to introduce Euler potentials globally for the magnetic field with 
non-trivial helicity, how to find the function ip, and why it has to be multi- valued, are not clear so 
far. For example, it is stated (Biskamp, 1993) that Euler potentials can not be introduced globally 
for a magnetic field with nonzero helicity unless the system is multiply connected. In (Sagdeev et 
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al., 1986) it is pointed out that magnetic field lines determined by the Lagrangian invariants do not 
admit any linkage, i.e., the helicity has to vanish. The representation ( |1.5| ) is used sometimes (Priest 
and Forbes, 1999; Wong, 2000) quite generally, but it is not mentioned that helicity has to vanish in 
this case, hence the structure of the magnetic field has to be relatively simple. 

The aim of this paper is to show how one can practically introduce Euler potentials ( |1.4| ) as well 
as the Clebsch representation ( |1.6| ) for a magnetic field with nonzero helicity. There is a solution 
of the MHD equations (Kamchatnov, 1982, Sagdeev et al., 1986) based on the Hopf invariant of 
the mapping of a 3D sphere into a 2D sphere S"^ (see, for example, Dubrovin et al., 1979). In 
this solution the magnetic flux tubes can link each other as many times as one wants. The MHD 
soliton has a known helicity following from topology, hence in each step of the calculation, there is 
opportunity to control the situation. Besides, this solution is relatively simple, and all the results 
can be obtained analytically. We will not use topological methods, because all our results can be 
obtained straightforwardly if some topological information has been taken into account from the very 
beginning. 

This paper is organized as follows. In Sections 2 and 3 we recall the details of the MHD 
Kamchatnov-Hopf solution. Euler potentials are introduced in Section 4. A visualization of the 
magnetic field structure is presented in Section 5. Helicity coordinates are introduced in Section 6, 
and Section 7 is devoted to the summary and discussion. 



2 MHD Kamchatnov-Hopf soliton 

First of all we will recall (Chandrasekhar, 1961; Kamchatnov, 1982) that any solinoidal vector field, 
divB = 0, gives rise to a solution of the steady-state MHD equations 



p(v.V)v = -VP+-^(B-V)B, (2.1) 

(v ■ V)B = (B ■ V)v, (2.2) 

divv = 0, (2.3) 

divB = 0, (2.4) 

in an incompressible plasma where the density p = const. Here P is the total (gas + magnetic) 
pressure, v is the plasma velocity. If we choose v = B/-^47rp, and P = const, then equations (pTII - 

O] ) are satisfied automatically. In this solution the magnetic tension is balanced by the centrifugal 



force. 

The idea of the Kamchatnov-Hopf soliton solution is to obtain a solenoidal vector field with known 
linkage using topological methods. A 3D sphere is defined in i?^ as a set of points {qi,q2, (13,(14) 
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such that qf + q2 + + qI = 1- Let us introduce two complex numbers Zi = qi + iq2, Z2 = qs + iq^, 
then can be described also as + = 1. A curve (a circle) on can be presented as 



\{t) = {Z^e'^^\Z2e 



iuJ2t'\ 



(2.5) 



where t is a parameter along the curve. It can be shown (Dubrovin et al., 1979) that two curves 
corresponding different initial points Zi, Z2 with integer numbers co'i,ci;2 link each other LJ1LJ2 times. 
A tangential field Y on 5*^ generated by the curve (|2.5|) is 



Y{uJi,UJ2) 



dl{t) 

dt 



-ujiq2,ujiqi, ~uj2q4,uJ2q3), 



(2.6) 



which also has the linkage ujiuj2- Now we can map the curve ( |2.5| ) into using the stereographic 
projection 



Xi 



qi 



1 + ^4 



, ^ = 1,2,3, 



(2.7) 



54 



qi 



z= 1,2,3. 



1 + X2' 

To obtain the vector field (|2.6|) in we can just differentiate equation 
parameter t 



(2.8) 

with respect to 



-{UJ2X1X3 + UJ1X2), {uiXi - UJ2X2X3), -UJ2{xl + xl - xl - 1) 



(2.9) 



Stereographic projection conserves the topological invariant that is the linkage ujiuj2. 

As a matter of fact, divJ 7^ 0, but it can be easily verified that the field B = 4 * J/(l + x^)^ is 
solenoidal, where x"^ = x\ + x\ + x\. The factor 1/(1 + x^)^ 7^ everywhere in i?^, therefore the field 
obtained. 



^ = TT"^ — r-:r{-2(u;22;iX3 + u;iX2), 2{ujiXi- UJ2X2X3), uJ2{xl + xl - xl - I)} , (2.10) 

(1 + x'^y^ 

has the same topological property as the field (|2.6| ) on S"^ . The factor 4 was introduced for the 
calculations convenience. 

The field ( ^.10] ) is the basis for the topological soliton. As was pointed out previously, if we 
introduce the plasma velocity v = B/v^47rp, and the pressure P = const, then MHD equations ( ^.1| - 
OD are satisfied automatically. We will refer to this solution as the MHD Kamchatnov-Hopf soliton. 
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3 Magnetic field lines. 



Let us now derive the equation of the magnetic field fines in R^. To tfiis end we can solve differential 
equations ^ = B using (|2.1CI|) , but it is mucfi more easy just to map tfie known integral curves (|2.5| ) 
from S*^ to witfi tfie fielp of stereografic projection (Kamcfiatnov, 1982; Sagdeev et al., 1986) 



Xl{t) 
X2{t) 
X3{t) 



2(xio cos(c<Jit) — X20 sin(co'it)) 
1 + Xq + (1 - Xq) cos(co'2t) + 2x30 sin(co'2t) ' 

2(x2o cos(c<Jit) + xio sin(ci;it)) 
1 + Xq + (1 - Xq) cos{uj2t) + 2X30 sm{iU2t) ' 

2x30 cos{uj2t) — (1 — Xq) sin(a;2t) 
1 + Xq + (1 - Xq) cos{uj2t) + 2x30 sin(u;2t) ' 



(3.1) 



Using trigonometric identities, it is possible to reduce equations p.l|) to tfie following form 



cos 61 

a + cos c)2 
sinOi 

X2 = —T 7^, 3.2 

a + cos c)2 



6 sin O 



3^3 



2 



a + b cos 62 
wfiere 



<di = Uit + ai, 62 = — (^2^ + ^2, CL = ~^ 



^ V 3-10 "I" 3-20 

1,2 2 1 4-^30 ~^ ~ -^o) -^10 

b = a — 1 = — f7^-2 2~^' cosai — 



20 



3:20 1 ~ 3;o . 2x30 o\ 

smai = — =, cos 02 = — , sma2 = — , (o.o) 



V3;?o + 3;io' V4xio + (l-3;2)2' V4xio + (1 

It turns out tfiat tfie magnetic field lines lie on tfie surface of tfie torus 



X, 



Xi = (a + 6 cos 62) cos 9i , 

X2 = (a + &COS 62) sin 9i, (3.4) 
X3 = 6 sin 02, 

wfiicfi is produced by tfie rotation of tfie circle x\ + (xi — a)^ = a? — 1 around tfie X3 axis. Tfie 
central torus degenerates into a circle (Sagdeev et. al, 1986), 

X3 = 0, x\ + x\ = 1, (3.5) 
wfiicfi will play an important role fiereafter. 
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4 Euler potentials 



It is convenient to choose as Euler potentials the following constants of integration (first integrals) 

of (OD 



{2ar 

or, in Cartesian coordinates. 



(4.1) 



X2 

a = UJ2 arctan \- ui arctan -, 

1 — 



Xi 



(1+X2)2- 

Then we can find the gradients of these functions 



(4.2) 



Va 



AuJiXiXs 



(x^-lY + xi 



UJ2X2 
2 i 2 ' 

1 2 



4OJ1X2X3 



(X^ 



1)2 + 



+ 



^0-2X1 
2 i 2 ' 



2co'i(xf + x| — X3 — r 



1)2 + x2 



(4.3) 



V/3 = - ^^^^,^3 [xi(x^ + x^-x^-l), X2(x2 + x^-x^-l), 2X3(X? + X^)], (4.4) 

and verify that equation (|1.4| ) is satisfied, i.e., the a,P are indeed Euler potentials. The potential 
a is a naked angle (i.e., an angle being not hidden under any trigonometric functions), which can 
have a nonzero contribution after integration of its gradient along a closed contour. Therefore, it is 
not surprising that first of all, a is a multi-valued function, and secondly, Va has a singularity on 
the circle (p.5|). 

The next step is to obtain the vector potential A = —pVa 



_ -4lJiXiX3(x^ +X^) UJ2X2 



R ' (X2 + 1)2' 

-AuJiX2Xz{xl + X2) CJ2X1 



R (x2 + 1) 

2uji{x\ + X2)(x^ + X2 — X3 — 1^ 
R 



(4.5) 



where R = (x^ + l)^((x^ — 1)^ + 4x3). Remember that the formal representation ( |1.5|) leads to 
zero helicity, but we know that K ^ hj the topological construction, hence, the potential ( |4.5| ) 
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should have some principal disadvantage. If vector potential A is defined by the differential equation 
(|1.2|), then we have to conclude that A is indeed the vector potential of the magnetic field B because 
the equation ( p. .21 ) is satisfied. But besides the differential equation there is also an integral equation 



(f{A-dl) = [ (B-dS) = Fb, (4.6) 
Jl Js 

where Fb is the magnetic flux through the surface S bounded by the contour L. Differential and 
integral equations sometimes are not identical, and in our situation this is exactly the case. If we 
choose any contour L which does not cross the disc bounded by the singular circle ( |3.5| ), then the 
circulation of A along L gives exactly the magnetic flux Fb- However, if the contour encounters the 
disc bounded by singular circle ( p.5|) , then the circulation gets an additional contribution 



jf(A-rfl) = FB + ^, (4.7) 

therefore, the integral equation (|4.6| ) is not satisfled. The formal reason for the multi-valued 
character of the circulation ( ^4.7|) lies in the singular behaviour of the latter at the circle ( |3.5| ), or due 
to the fact that the function a ( |4.1| ) is a naked angle. Hence, we have to proceed with the Clebsch 
representation (|1.6| ), and to flnd a function ip to compensate the singularity in the potential ( |1.5| ). It 
is clear that the function ip has also to be a naked angle like the function a, and its gradient should 
have a singularity at the circle ( ^3|) 



1 — 1 1 71 

ip = -ui arctan(— ) = -uji{-uj2t + ^2 + tt)- (4- 

4 2x3 4 2 

Then, the Clebsch potential (p.. 61) turns out to be 



r U02X2 + 1^1X1X3 (^1X2X3 - U02X1 (1 + X 



~\ (1 + X2)2 ' (1 + X2)2 ' 2(1 + X2)2 J' ^'''^ 

It has no singularity in the whole space like the magnetic fleld ( 2.101 ), and both the differential 
(|1.2|) and the integral ( [4.6|) equations are now satisfled. 

It is interesting to note that the Clebsch representation ( |1.6|) formally looks similar to the gauge 
condition A' ^ A + V0. Nevertheless there is an essential difference. The function has to be a 
single- valued one for the gauge transformation at least for the simple connected region Q, hence the 
integral of its gradient along any closed contour has to vanish. Contrary, the function ip in the Clebsch 



representation ( |1.6|) has to be a multi-valued one, and the integral of its gradient along some closed 
contour can have nonzero contribution. Generally speaking, the question whether the gauge function 
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is a multi or single valued one is not really important for many applications in electrodynamics. But 
for such a delicate characteristics of the field as the magnetic helicity, the solution of this question 
plays the key role. It is the multi-valued function ip which does the nonzero helicity. 

Using the vector potential (|4.9| ) and the magnetic field (|2.1CI|) , we can calculate the helicity as the 
volume integral (|L1|) 



/<- = -^. (4.10) 



The negative sign in ( |4.1CI| ) is connected with the parameter t in the initial curve ( p. 51 ) at , so 



that e*^* gives a minus, whereas e gives a plus. 

The Clebsch representation ( |1.6|) leads to another way to compute the helicity 



ir = W /(B- dS), (4.11) 



s 



which, of course gives the same result ( |4.1(]|) . Here S is the singular circle (|3.5| ), and {ip} is the 



jump of the function ip on the latter. As one can see, helicity can be calculated from the surface 
integral ( ^.111 ) rather than from the volume integral ( |1 . 1| ) , which is simpler to do. It is also interesting 
that the helicity is equal to the magnetic flux through the singular circle times the jump of the ip 
function. 



5 Visualization 

It is worthwhile to present pictures of the magnetic field structure of the MHD Kamchatnov-Hopf 
soliton solution as mathematical examples for illustration. 

We start with the simplest case, cui = 0^2 = 1 . The flux tube looks like a torus twisted by the 
angle 360°. To see this more clearly, the tube presented is chosen to have a rectangular cross section 
(Figure 1) so that one can easily follow the screwed color boundaries. 

The surface Euler potential [3 = const is just a usual torus (Figure 2), it stays more or less the 
same for all uji,uj2- The magnetic field lines are swept around this torus. 

The surface a = const is more complicated (Figure 3). It is similar to a ribbon twisted by 360". 
Such a surface can not be continued to the closed one in without self crossing which is because 
Va has a singularity at the circle (|3.5|) . 

There is a simple way to imagine the magnetic field structure. Let us take a paper ribbon, twist 
it by the angle 360° (note that twisting by 180° gives a Moebius sheet), glue the edges of the ribbon 
together, and then cut it along the central line with scissors. As a result, we get two ribbons linked 
to each other. If we continue this procedure and cut the two ribbons obtained along their central 
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axis, and so on, we can observe that each ribbon hnks any other one exactly one time (Figure 4). 
This behaviour is reflected in the topological invariant helicity K ( ^.101 ). 

It is difficult to imagine that the intersection of two surfaces a = const and (3 = const for different 
constants can give linked lines, nevertheless, it is so. 

To complete the case cji = = 1, we present also the surface ip = const (Figure 5) which has a 
spiral structure converging to the singular circle (|3.5|). 

After these relatively simple pictures we can proceed to the general case. First we recall that 
two numbers are relatively prime, if and only if the greatest common divisor of the numbers, is 
one. For integers uji,uj2 such that Ui = n,uj2 = m are relative prime, the magnetic field lines of the 
MHD Kamchatnov-Hopf soliton are linked into (n, m) knots which are topologically nonequivalent 
for different {n,m). They form the known family of toric nodes (Crowell and Fox, 1963). 

The case cji = 2,u;2 = 1 is depicted in Figure 6 (single flux tube). Figure 7 (surface a = const), 
and Figure 8 (knot (2,1) ). The more complicated case tui = 2, u;2 = 3 is presented in Figure 9 (single 
flux tube). Figure 10 (surface a = const), Figure 11 (central fragment of the surface a = const), and 
Figure 12 (knot (2,3) ). 

It is interesting that the surface of the Euler potential a = const for the latter case (Figures 10, 
11) is similar to a propeller, and this circumstance seems not to be a pure coincidence. The propeller 
has to create curls of air for producing a moving force, and at least some surfaces a = const (Figures 
3, 7, 10) might be used for this aim just from topological reasons. Of course there is the question 
about the efficiency of such airscrews or waterscrews, but this is not a subject of this paper. 



6 Helicity coordinates 

A magnetic field line is defined by two Euler potentials a, (3, and a point on this line is controlled by 
the parameter t. We can use another parameter along the magnetic field line instead of t. Then 
a,P,ip, i.e., all functions taking part in the Clebsch representation of the vector potential ( |1.6|) , can 
be used as new curvilinear coordinates which have some useful property. 

We already have an expression for the Clebsch coordinates via Cartesian coordinates ( |4.1| , [4.2| , 
Of ). It is possible to simplify these equations noting that without loss of generality, we can assume 



ai = in (11) and then obtain 



2X3 

a = oji arctan 



l-a;2' 



„2 



, 1 ,x'-l. 
ip = -uji arctan( 



4 ' 2x 



3 
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Now we can also find the mapping (xi, X2, x^) (a, P, ip) 



Xi 



v/r^cos(M-|) + i' 



2v^sin(^ + ^ - ^) 

X2 = J -^ (6.2) 

v/r^cos(M-f) + l' 

v/T^sin(M - I) 



2:3 



v/r^cos(S^-f) + i 



1^1 



After some algebra one can find the Jacobian of this transformation 



_ D{Xi,X2,X3) _ A 



where 



A = 8(v/(l-4/?) + 4/^7 V(l - 4/?) + 2(1 - 4/^)572 - 8/^7 + 27 - ^(1 - 4/?)72), 
5 = uj2UJi{--f^ + 12/^7^ - 48/5^75 + 64/?=^7^ - 5(1 - - 107^ + 80/57^ - 160/3 V 

- 10(1 - 4/3)572 - 57 + 20/37 - V(l - 4/3)), (6.4) 

7 = cos(-4— + -). 

This equation is a bit complicated, nevertheless it is possible to verify that J 7^ in the whole 
space, hence the coordinates (a,/3, ■?/') can be introduced in R^. 

Let us compute the magnetic helicity using these new coordinates 



K = (A ■ B)dxi A dx2 Adx3 = / {Vip ■ [Va x VP])dxi A dx2 A dxs = 

dxi A dx2 A dx^ = / dadpdijj, (6.5) 



^3 D(a;i,a;2,a;3) 



where it is supposed that the whole space is mapped into the region Q', —>■ Q' under the 
transformation {xi,X2,X3) — > {a,(3,ip). Therefore, it turns out that in the new variables, magnetic 
helicity is equal to the volume of the configuration space (a, /3, ■?/'). It is easy to verify that the new 
formula ( |6.5| ) gives the same result ( [4.10 ) for the helicity if we take into account that the coordinates 



[a, /?, ip) are varied within the following limits 



ae(-7r,7r]; /3g[0,^); ip e {-^uJiUJ2TT,^uJiUJ2rc], (6.6) 
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where brackets ( or [ are used to show that the element close to the bracket is excluded or included 
in the list of elements, respectively. 

One can see that the space is mapped onto the parallelepiped ( |6.6D in which the straight lines 
[a = const, (3 = const) represent the magnetic field lines. It is surprising that all the complicated 
magnetic structure is converted into a very simple geometrical object, that is the parallelepiped (|6.6|) . 
In fact, the situation is not that simple. To make field lines which are closed, we have to glue the 
end points. The points on the left boundary a = — tt have to be considered identical with those on 
the right boundary a = n, and the points of the bottom boundary ip = —jUiU2iT are identical with 
those on the upper boundary ip = jU!iU!2T!' after the rotation of the latter by the angle 2c(JiC(j2vr. 

We note that helicity coordinates can be particularly important for numerical simulations where 
helicity control is required. 



7 Discussion and summary 

It was shown that Euler potentials can be introduced globally for a magnetic field with nonzero 
helicity even for the simply-connected region (the space in our case), contrary to the remark of 
Biskamp (1993). Therefore most of the coordinate systems (Pudovkin and Semenov, 1985; Pustovi- 
tov, 1999) based on Euler potentials ( p..4|) , such as the helicity system ( |6.2|) can still be applied also 
to magnetic structures with nonzero helicity K ^ 0. On the other hand, one has to be particularly 



careful with the vector potential. Remember that the simple representation ( |1.5| ) is not appropriate 
for the magnetic field with K 0, instead the Clebsch representation (|1.6|) has to be used. 

As we saw, the function plays a key role in calculating the magnetic helicity. The multi-valued 
character of this function is connected with the singular behaviour of the gradient of at least one 
Euler potential {a in our case). In its turn, the singularity of the Euler potential is the consequence 
of the fact that the a = const surface is highly twisted for the case K and cannot be continued 
to a closed surface in R^. 

The helicity turns out to be equal to the magnetic fiux through the singular circle times the jump 
of the function ip, hence, the calculation of K can be reduced to a surface integral.lt seems that the 
simple formula ( [4.11| ) can be extended to the general case as 



K = J2{^AFbj, (7.1) 
j 

where Fbj is the magnetic fiux through the surface bounded by j-singular curve of the Euler 
potential and {ipj} is the jump of the function ip at this surface. 

If by chance it is known that all singular lines of the Euler potentials lie on a surface S (the 
surface of the Sun, for example) then the magnetic helicity can be found using only data of the 
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normal component of the magnetic field -B„ on S from the equations ( 4.11 or which is an 
important problem for solar physics. But if a singular line is inside the Sun, it is not possible to find 
the helicity using surface data. The maximum of what can be done is to estimate the helicity if one 
could somehow control the magnetic flux closed under the Sun's surface. 

The Kamchatnov-Hopf solution seems to be the simplest one which can describe the magnetic 
field with such a non-trivial helicity. Therefore it may play the same role for the investigation of 
different helicity problems as the Harris (1963) layer in plasma physics or the Petschek(1964) solution 
in reconnection theory. 



12 



Acknowledgements: We are grateful to I. V. Kubyshkin for helpful discussions. Part of this work 
was done while VSS and DBK were on a research visit to Graz. This work is partially supported by 
the Russian Foundation of Basic Research , grant No. 01-05-64954, by the INTAS-ESA grant No. 
99-01277, and by the programme INTERGEOPHYSICS from the Russian Ministry of Higher Edu- 
cation. Part of this work is supported by the Austrian "Ponds zur Forderung der wissenschaftlichen 
Forschung" , project P13804-TPH. We acknowledge support by the Austrian Academy of Sciences, 
"Verwaltungsstelle fiir Auslandsbeziehungen" . 



13 



References 

[1] Biskamp, D., Nonlinear magnetohydrodynamics, Cambridge University Press, 1993. 

[2] Chandrasekhar, S., Hydrodynamic and hydromagnetic stability, Oxford University Press, 1961. 

[3] Crowell, R. and R. Fox, Introduction to knot theory, New York, 1963. 

[4] Dubrovin, B. A., S. P. Novikov, and A. T. Fomenko, Modern geometry, Nauka, Moscow, 1979. 

[5] Harris, E. G., On a plasma sheath separating regions of oppositely directed magnetic fields, 
Nuovo Cimento, 23, 115, 1962. 

[6] Kamchatnov, A. M., Topological soliton in magnetohydrodynamics, Sov. JETP, 82, No 1, 117, 
1982. 

[7] Moffat, H. K. , Magnetic field generation in electrically conducting fluids, Cambridge University 
Press, 1978. 

[8] Petschek, H. E., Magnetic field annihilation, NASA Spec. Publ, SP-50, 425, 1964. 

[9] Pudovkin, M. I. and V. S. Semenov, Magnetic field reconnection theory and the solar wind- 
magnetosphere interaction: A review. Space Sci. Revs., 41, 1, (1985). 

[10] Priest, E. and T. Forbes, Magnetic reconnection, Cambridge University Press, 2000. 

[11] Pustovitov, V. D., Magnetic coordinates with double straightening. Plasma Phys. Rep., 25, No 
12, 963, 1999. 

[12] Sagdeev, R. Z., S. S. Moiseev, A. V. Tur, and V. V. Yanovskii, Problems of the theory of 
strong turbulence and topological soliton, in Nonlinear Phenomena in Plasma Physics and 
Hydrodynamics, edited by R. Z. Sagdeev, (MIR Publishes, Moscow, 1986), p. 135. 

[13] Schindler, K., M. Hesse, and J. Birn, General magnetic reconnection, parallel electric fields and 
hehcity, J. Geophys. Res., 93, No A6, 5547, 1988. 

[14] Taylor, J. B., Relaxation revisited, Phys. Plasmas, 7, No 5, 1623, 2000. 

[15] Wiegelmann, T. and J. Biichner, Evolution of magnetic helicity in the course of kinetic magnetic 
reconnection. Nonlinear Processes in Geophysics, 8, No 3, 1623, 2001. 

[16] Wong, H. v.. Particle canonical variables and guiding center Hamiltonian up to second order in 
the Larmor radius, Phys. Plasmas, 7, No 1, 73, 2000. 



14 



[17] Yoshizawa, A., N. Yokoi, and H. Kato, CoUimation mechanism of magnetohydrodynamic jets 
based on helicity and cross-helicity dynamos, with reference to astronomical jets, 

Phys. Plasmas, 7, No 6, 2646, 2000. 



15 



Figure Captions 



1. 


Magnetic flux tube for the case a;i = 1, a;2 = 1. 




2. 


Surface (3 — const for the case a;i = 1, a;2 = 1. 




3. 


Surface a — const for the case uji — 1,uj2 — 1. 




4. 


Two hnked flux tubes for the case uji = l,ui2 = 


1. 


5. 


Surface ■0 = const for the case cui — 1,002 — i- 




6. 


Magnetic flux tube for the case cui — 2,u!2 — 1- 




7. 


Surface a — const for the case cui = 2,002 ~ 1. 




8. 


Two hnked flux tubes for the case loi — 2,u;2 — 


1. 


9. 


Magnetic flux tube for the case a;i = 2, a;2 = 3. 




10. 


Surface a — const for the case a;i = 2, a;2 = 3. 




11. 


Central fragment of the surface a — const for the case cui — 2,lo2 — 3 


12. 


Two hnked flux tubes for the case a;i = 2, a;2 = 


3. 
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